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Abstract. Combining XMM-Newton and Chandra data, we have performed a detailed study of Abell 3558. Our analysis shows 
that its dynamical history is more complicated than previously thought. We have found some traits typical of cool core clusters 
(surface brightness peaked at the center, peaked metal abundance profile) and others that are more common in merging clusters, 
like deviations from spherical symmetry in the thermodynamic quantities of the ICM. This last result has been achieved with 
a new technique for deriving temperature maps from images. We have also detected a cold front and, with the combined use 
of XMM-Newton and Chandra, we have characterized its properties, such as the speed and the metal abundance profile across 
the edge. This cold front is probably due to the sloshing of the core, induced by the perturbation of the gravitational potential 
associated with a past merger. The hydrodynamic processes related to this perturbation have presumably produced a tail of 
lower entropy, higher pressure and metal rich ICM, which extends behind the cold front for 500 kpc. The unique characteristics 
of A3558 are probably due to the very peculiar environment in which it is located: the core of the Shapley supercluster. 
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1. Introduction 

Superclusters are the most massive objects in the Universe and 
are destined to collapse. By studying superclusters, we can ob- 
serve the processes related to individual components at an early 
stage of merging. These are: mutual interaction between the 
components, accretion history of the member clusters, interac- 
tion between the collapsed objects and the intra supercluster 

medium. 

The Shapley Supercluster ( Shaplevl[l930l) is one of the largest 



concentrations of mass i n the local Unive rse dBardelli et al 
19941; IZucca et all 1 19931; lEttori et all 1 19971). The core of the 



concentration is dominated by the greatest cluster complex of 
the supercluster, which is composed of A3562, SC 1329-313, 
SC 1327-312, A3558 and A3556 (from SE to W). Many mul- 
tiwavelength observations have been carried out on the core of 
the Shapley Supercluster, from which the total gravitating mass 



pansion and contract, even if smaller than that needed to virial- 
ize. Therefore, they suggest that the whole A3558 complex is 
now in an early phase of contraction. 

The central galaxy cluster in the core region is A3558, which is 
also the brightest and most massive cl uster in the complex. The 
ROSAT image of the core complex dKull & Bohringeii Il999h 
shows that the X-ray surface brightness of A3558 is elongated 
in the SE-NW direction, along which th e other members an d 
the galaxy distribution are also aligned dBardelli et al.lll994l) . 
Therefore the infall of matter onto A3558 should be dominant 
along this direction. 

The temperature distributio n of the ICM has been investi - 
gated using ASCA data bv iMarkevitc h & Vikhliiiinl dl997h . 



Hanami et alj (Il999h and lAkimoto et al., (2003») . showing de- 
viations from spherical symmetry. 

A radio analysis of the A3558 complex (IVenturi et al. I I2OOOI) 



is estimated to be 10 



10 



16 1,-1 



(iMetcalfeet al.Ml994 



Ettori et al.lll997h . The detected low surface brightness X-ray 
emission connects A3558 wit h the two groups and A 3562 
and probably also with A3556 dKull & Bohringeilll999l) . Also 
the galaxy iso density contours extend over the whole complex 
dBardelh et al.., 1998.), indicating th e prese nce of a structure on 
supercluster scales. lAkimoto et al.l (l2003b have calculated the 
mean density of the core region of the Shapley supercluster 
to be 25 times greater than the critical density. This density is 
larger than the value needed to decouple from the cosmic ex- 



underlines the absence of extended radio sources (halos and 
reUcs, which are observed in many merging clusters) in A3558. 
Also the radio-optical luminosity function for early-type galax- 
ies in the complex is significantly lower than that deiived for a 
sample of early-type galaxies in clusters. 
In this paper we present results from XMM-Newton and 
Chandra observations of A3558. In Section 2 we present the 
analysis and describe a new technique to produce temperature 
maps; in section 3 we present the most interesting features of 
the cluster and in section 4 we summarize our results and dis- 
cuss the dynamical state of the cluster. 
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Table 1. Source and background observations. 



2.2. Chandra data reduction 



Obj name 


Nh 


SXR bkg,;^ 




(1022cm-2) 


(cts s^'arcmin"^) 


A3558 


0.0389 


(301 ± 10)10-^ 


PHL 5200 


0.0526 


(111 ± 8)10-" 


PKS 0312-770 


0.0797 


(122 ±6)10-"^ 



" Average X-ray background count rate in the 3/4 keV band, derived 
from the ROSAT all-sky survey diffuse background maps. 



At the nominal redshift of A3558 (z — 0.047), 1 arcmin corre- 
sponds to 54 kpc (//() = 70/170 km s"' Mpc"'). All quoted errors 
are at lo", unless otherwise stated. 



2. Data analysis 

2.1. XMM-Newton data reduction 

A3558 was observed by XMM-Newton during revolution 388 
with the EPIC detectors in Full Frame Mode with the thick 
filter We have obtained calibrated event files for the MOSl, 
MOS2 and PN cameras with SAS version 6.0. These files 
are then manually screened to remove remaining bright pixels 
or columns. To reject soft proton flares, we accumulate light 
curves in the 10-12 keV energy band, where the emission is 
dominated by the particle-induced background. A visual in- 
spection of the light curve reveals that the observation is not 
badly affected by soft pr otons: we remove only 0.5 ks, using 
the cr-clipping technique dMartv et al.ll2002h . leaving an expo- 
sure of 43 ks for the MOS cameras and 35 ks for the PN. 
For the study of galaxy clusters, especially in their outskirts, it 
is very important to subtract the instrumental background. The 
quiescent particle induced background can be properly sub- 
tracted using a large collection of background data, such as 
long observations of blank sky fields. In this case, the choice 
was not easy because there are not many blank sky fields ob- 
servations performed with the thick filter We have selected two 
fields PHL 5200 & PKS 0312-770 which have similar charac- 
teristics to the source (see Table [TJ. We have applied to these 
observations the same selection criteria applied to the source, 
we removed point sources, reprojected the files to the observa- 
tion coordinates and merged the fields with the SAS task merge, 
for a total exposure of 79 ks for the PN and 8 1 ks for the MOS 
detectors. 

The intensity of the quiescent particles background may differ 
during the observations of the source and of the blank fields. 
We have compared the count rates at energies greater than 10 
keV, where the emission is dominated by the background, in 
the source and blank sky fields: the level of count rates during 
background observations is about 10% higher than that of the 
source observation. To correct for this effect we renormalized 
the background by multiplying it by a factor of 0.87 before per- 
forming background subtraction. 



Chandra observed A3558 in 2001 with the Advanced CCD 
Imaging Spectrometer (ACIS). We have created standard level 
two event files, screened for background flares using version 
2.4 of the ACIS calibration products available with the version 
3.02 of the CIAO package. The useful time for the analysis af- 
ter the screening is 13.9 ks. 

We have performed background subtraction using the "blank 
sky" datasets enclosed in the CIAO Calibration Database. We 
have processed the background file with the same procedure 
of the source file and we have reprojected it to the Abell 3558 
observation sky coordinates. 

2.3. Maps from images 

Using data from X-ray imaging spectrometers, such as EPIC 
and ACIS, we can derive a 2-dimensional temperature map of 
the cluster. In order to exploit the spatial resolution of these in- 
struments, we want to obtain at the same time as much "spec- 
tral" and "spatial" information as possible. To do this, we have 
derived temperature maps from X-ray images using two dif- 
ferent and comple mentary te chniques: the wavelets -H hard- 
ness ratio analysis ( Briel et al.li2004) and the broad band fitting 
with adaptive binning, that we will describe in section 12.3.11 
Combining the temperature map and the surface brightness 
image, we also derived pseudo-pressure and pseudo-entropy 
maps. Besides giving informations themselves, the maps can 
be used to select interesting regions for a proper spectral anal- 
ysis. 

2.3.1 . Maps from adaptive binning + broad band fitting 

Temperature maps of galaxy clusters are o ften obtained ap 



plying adaptive smoothing to X-ray images (iMarkevitch et al 
■2000) . in order to emphasize spatial details. The smoothing 
technique increases the local signal to noise ratio by correlat- 
ing neighboring pixels through a convolution process. One of 
the problems of this technique is that it is often difficult to han- 
dle the statistics of the correlations introduced by the smooth- 
ing. The adaptive binning technique locally groups and aver- 
ages together data, allowing us to have a firmer control on the 
statistics: this is why it should be preferred when a precise error 
estimation is needed. 

Rather than using the hardness ratio technique, we have per- 
formed a "broad band fitting" starting from images, as in 
Markevitch et al. I (l2000h . The combined use of these techniques 



allows us to have a good control on the statistics and to produce 
quantitative thermodynamic maps (i.e. maps with their associ- 
ated errors). In the following, we will describe the application 
of this technique to EP IC images, but it can also be easily ap- 
plied to ACIS data (see lGhizzardi et al ]l2005l for an application 
to Chandra images). 

We extract images in five energy bands (0.4-0.8 keV, 0.8-1.2 
keV, 1.2-2 keV, 2-4 keV, 4-10 keV) from the event files of the 
observation and of the background and we generate exposure 
maps in the same energy bands with the SAS task eexpmap. 
Then we select the images with lowest signal (i.e. one of the 




Fig. 1. Maps of the main thermodynamic quantities of the cluster, obtained with the adaptive binning and broad band fitting 
technique. Upper left: "projected emission measure" (XSPEC NormaUzation per pixel); Upper right: temperature (keV), Lower 
left: pseudo-pressure (arbitrary units); Lower right: pseudo entropy (arbitrary units). Coordinates on the maps are right ascension 
and declination. 



MOS images in the hardest energy band) and we perform adap- 
tive binning on this image, in order to group pixels in bins with 
the same signal to noise ratio. To do this we have used the 2- 
dimensional adaptive binning technique bv lCappellari & Copin 
(120031) . which is based on the Voronoi tessellation. This pro- 
cedure has been applied to the MOSl image of A3558 in the 
4-10 keV energy band with a target signal to noise ratio of 10; 
the same binning is applied to all images, assuring a minimum 
S/N=10 in all bins in all images, since the statistics in the PN 
and in the other bands is better than in the image we have se- 
lected. For each bin, we sum the count rates of MOSl, MOS2 
and PN, obtaining the EPIC count rate and its error in each band 
and we compare them with a grid of predicted count rates, gen- 
erated with XSPEC using absorbed mekal models with differ- 
ent temperatures and normalizations, while the absorption col- 
umn and redshift are fixed to the known values for this cluster 
and the iron abundance is fixed at 0.3 solar We select our "best 
fit" values for temperature and normalization as those that min- 



imize j^'. Using the;t' matrix, we can also find values within a 
fixed Ax^ and we can associate errors with best fit parameters. 
The normalization values are in XSPEC units: 



K = 



10 



14 r- 



rindV, 



(1) 



where Da is the angular size distance to the source (cm), 
is the electron density (cm"-'), and «h is the hydrogen density 
(cm^^). If we divide the normalization by the number of pixels 
in each bin {K/Np), we obtain the "projected emission mea- 
sure" that is proportional to the square of the electron density 
Hp, integrated along the line of sight. 

Using this "projected emission measure" we can obtain for 
each bin the "pseudo-entropy" s and "pseudo-pressure" P, de- 
fined as: 



Nn 



-1/3 



P = T\ — 



1/2 



(2) 
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Fig. 2. Center of A3558 (inner 160 arcsec). Upper panel: 
Temperature map, lower panel: map (3 d.o.f.). Coordinates 
on the maps are right ascension and dechnation. 



and we can produce maps of these important thermodynamic 
quantities. 

The maps obtained with this technique with EPIC are shown 
in Fig. [T] Our technique can also produce for each of these 
maps, two maps corresponding to the upper and lower limit at 
a given confidence level. These maps are extremely useful to 
determine whether the gradients that we observe in Fig. [T] are 
significant or not. We have applied this analysis to the "blue 
spot" in the temperature map at the center of the cluster (Fig. 
|2^), corresponding to the brightest bin. Its temperature, kT= 
(3-7tQ jg) keV, is significantly lower than the mean value of the 



adjacent bins, kT=(4.89 + 0.11) keV. Moreover the of the 
"fit" in this bin is higher than in the surrounding bins (Fig.|2]3), 
suggesting that the fit in the blue spot is worse than in the other 
regions. Investigation of the Chandra image (Figl?) shows that 
the blue spot coincides with the AGN located in the BCG. This 
shows that our technique is able to recognize the spectral shape, 
even if we use large energy bands. 

The typical relative errors on the temperature are of the order 
of 3 - 5% in a circle of radius 7 arcminutes from the SB peak, 
while at outer radii, errors grow to about 5 - 10%. The larger 
errors (~ 10%) are found in the southern and western regions, 
where the surface brightness decreases more rapidly: in these 
regions the dominant factor of indetermination in temperature 
is the background subtraction. 

We do not show here the maps obtained from ACIS images: 
because of the short exposure of this observation we cannot 
produce maps with bins smaller than in EPIC maps. A deeper 
Chandra observation could allow us to exploit the full potential 
and the spatial resolution of the ACIS instrument. 

2.3.2. Maps from wavelet + hardness ratio 

In Fig. [3] we show the maps of the main thermodynamic quan- 
tities obtained with hardness ratio (HR) and wavelet decom- 
position, with the technique described in Briel et al. (2003). 
The interest in this technique is more focused on detecting the 
regions where there are temperature gradients than in the ab- 
solute values, since they are used to select regions for a proper 
spectral analysis. The main advantage of the wavelet decom- 
position is to reveal the small-scale features, which may not be 
resolved by the binning method. 

A qualitative comparison with the maps in Fig. [T] shows that 
they are consistent: the features detected with the new tech- 
nique are seen also in Fig.|3] 

Only the temperature maps can be compared quantitatively 
(Fig. [TJi and Fig. [3}5), since the other maps have different 
units. The fact that temperatures are always higher in the HR 
map is not worrying because the purpose of this technique is 
not to produce quantitative temperature maps. Moreover, the 
simplified HR technique that we have used here assumes that 
the observation is performed with a medium filter and therefore 
can overestimate the temperature. 

2.4. Maps from spectral analysis 

In order to test our technique we have compared the tempera- 
ture and normalization maps that we derive from images with 
the same maps obtained from a proper spectral analysis. With 
the results of the spectral analysis we have produced tempera- 
ture and "projected emission measure" maps that can be com- 
pared to the maps obtained from images (FiglT]). 
A plot of the temperature obtained with our technique versus 
the temperature obtained with spectral analysis is shown in Fig. 

a 

The mean absolute deviation {\Ti,„agelT spectral - 1|) is 8.8%, 
with a 1.1% dispersion around the mean. The largest deviation 
between the temperature estimates is (T/,,, - Tsp) = 4.5(0"^^^ + 




Fig. 3. Maps of the main thermodynamic quantities of the cluster, obtained with wavelet decomposition and hardness ratio. Upper 
left: wavelet decomposed image of the cluster; Upper right: temperature (keV), Lower left: pseudo-pressure (arbitrary units); 
Lower right: pseudo-entropy (arbitrary units). Coordinates on the maps are right ascension and declination. The hot zone in the 
NW is an artefact due to both high background in the hard band and to its position near the edge of the FoV. 



crjp^)"-^). A similar analysis of the difference between the tem- 2.5. Ratio maps 
perature obtained using MOS 1 data and that obtained from the 
PN shows comparable variations (mean deviation 17.3% and 
scatter 5.3%). 

Therefore the error associated with our technique is some- 
what smaller than the systematic error associated with cross- 
calibration uncertainties between MOS and PN. 
The good agreement between the results of our technique and 
those of the spectral analysis confirms that the maps obtained 
from images can substitute for those obtained from the spectra 
in all the cases where a proper spectral analysis would require 
too much time, as for instance in long observations of bright 
clusters or in large samples. 



The most interesting features in a non-relaxed galaxy cluster 
are the deviations from the spherical symmetry, which would 
be expected in an ideal relaxed cluster Such deviations can be 
seen in any ICM observable (see for instance the "residual" 
image of the Coma Cluster, by Neumann et alj|2003 ). We have 
decided to look for deviations in pressure (since it is strongly 
related to the Dark Matter distribution) and in entropy. First 
we have derived pressure and entropy in all the polygons pro- 
duced by the Voronoi tessellation, using the results of spec- 
tral analysis and m aking the volume es timates as described in 



Henry et al.l (|2004|) and iMahdavi et all (120051) . Then we cal- 



culate the predicted "relaxed" values, using the best fit mod- 
els t o the pressure and entr opy profiles in a sample of clus- 
ters (" Finoguenov et alJ2005 ) and centering them on the surface 
brightness peak (the impact of the choice of the center of sym- 
metry will be discussed later). Finally we divide the observed 
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Fig. 5. Ratio maps of pressure (left) and entropy (right). The regions labeled "*" have no data because of the subtraction of point 
sources. The white circle is the region considered for the analysis of the pressure excess, while the black circle is the reference 
region. Coordinates on the maps are right ascension and declination. 




4 6 
Spectral Temperature (keV) 



Fig. 4. Temperature derived in each polygon with our technique 
versus the weighted mean of the best fit temperature in each 
detector obtained from spectral analysis. The continuous line 
represents the ideal expected relation Tj^ag^, - T spectral- 



values by the predicted one in each polygon, obtaining the ratio 
maps shown in Fig. |5] The area-weighted fractional scatter in 
the entropy and pressure is 0.23 + 0.01 and 0.27 ± 0.01, respec- 
tively, inside a circular region of r = O.Srgoo (r ^ 9.2 arcmin 
from the cluster center). 

Since this technique is strongly based on geometrical and sym- 



metry assumptions, we have recalculated the maps using an- 
other point as the center of the cluster This "new center" is the 
center of the surface brightness isophotes (we will discuss in 
Sec.|3]the choice of this point) and is located ~ 70 arcsec SE of 
the surface brightness peak. In Table |2] we compare the frac- 
tional scatter obtained using the different centers. 
The choice of the center does not have a strong impact either 
on the values of the pressure and entropy scatter, especially at 
large scales, or on the ratio maps. The structures that appear in 
Fig. |5] are seen also in the maps with the other center (that we 
do not show here). Some obvious changes can be found only 
in the inner regions: if we use the "SB center" the main surface 
brightness peak appears as a substructure. 

3. Results 

The morphology of the cluster in X-rays can provide important 
information on its dynamical state. In Fig. |6] and |7] we show 
the images from EPIC and ACIS-S3 respectively. The peak of 
the surface brightness coincides with the BCG, which hosts an 
AGN. Comparing the imag e of this cluster with i mages from a 
sample of nearby clusters dGhizzardi et al.ll2005b we have ob- 
served that the X-ray emission is not as peaked as in cool core 
clusters, but the peak is more apparent than in the clusters that 
are currently undergoing a major merger. 
The surface brightness contours in Fig.|6]show that the distribu- 
tion of the ICM is elliptical, elongated in the SE-NW direction 
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Table 2. Fractional entropy and pressure scatter obtained using 
the surface brightness peak (SB peak) and the center of large 
scale X-ray isophotes (SB center) as center of symmetry. 



Region^ 


center 


Entropy 


Pressure 


r < O.lrsoo 


SB peak 


(0.179 ±0.011) 


(0.629 ± 0.011) 




SB center 


(0.193 ±0.008) 


(0.667 ± 0.012) 


r < Q.lrsQo 


SB peak 


(0.180 ±0.006) 


(0.413 ± 0.006) 




SB center 


(0.161 ±0.005) 


(0.448 ± 0.007) 






m 907 J- n nrm 


in 11 f, J- n r\nQ\ 






a.i 


1.0 


10.0 




Fig. 6. This EPIC flux image in the 0.4-2 keV band is com- 
puted following a procedure similar to the one described in 
Baldi et al.l (120021) . The EPIC source counts image, S epic, is 



computed by summing the MOSl, MOS2 and PN source im- 
ages, while an EPIC source exposure map, Iepic, is computed 
by summing the source exposure maps of each detector. The 
EPIC count rate image is defined as: crgpic = S epic /(epic - The 
total conversion factor (count rate to flux) c/epic has been cal- 
culated using the exposure times for MOS 1 , MOS2 and PN, the 
conversion factors for the three instruments, c/^osi, c/mosi, 
cfpN, and following the formula: = + + 
It is easy to show that the EPIC source flux image is then : 
Fepic - c/epic ■ ctepic- Similarly we compute the EPIC back- 
ground flux image from the MOSl, MOS2 and PN background 
images and exposure maps. Finally by subtracting the EPIC 
background flux image from the source flux image we derive 
an EPIC net flux image in units of 10^'^ergs cm"^ s^'pixel"' 
(one pixel is 8.5 * 8.5arcsec^) . Coordinates on the image are 
right ascension and declination. 



(inclination= 50°). Excluding the inner regions (r < 120 arcsec 
from the SB peak), the cluster appears remarkably symmetric 
large scales. Approximating the contours in Fig|6]with ellipses, 
we found that their center does not change significantly if we 
exclude the inner contours. This center can be considered the 
center of symmetry of the cluster: it is located ^ 70 arcsec SE 
of the surface brightness peak. Such offsets are not uncommon 
in clusters. 

In the Chandra image (Fig.|7]l there is a clear brightness discon- 
tinuity (visible even if at a lower resolution in Fig|6]l in the NW 
sector. This edge is similar to those observed with Chandra 
in many clu sters dM cAevitch et al.l200oH IVikhlinin et al ]|2001 



Dupke & W hite 2003|): the analysis of :he surface brightness 
profile and of the temperature discontinuity (Section [3. Il l con- 
firms that it is a cold front. 

The temperature maps (Fig. [TJ) and Fig. [3J)) confirm that the 
cluster is not spherically symmetric: we see a "hot bridge" 
(kT - (6.2 ± 0.5) keV) which surrounds the cluster center on 
three sides, while the temperature in the SE sector at the same 
distance from the center is significantly lower {kT - (5 . 1 ± 0.4) 
keV, ~ 20%), similar to the temperature of the cluster core. 
The entropy maps (Fig.[TJl and Fig.[3}l) do not show the spher- 
ical symmetry typical of relaxed clusters: the low entropy gas 
extends from the core towards the SE for about 500 kpc. The 
pressure maps (Fig. [1]; and Fig. |3};) seem to suggest that this 
quantity is more spherically symmetric than the others but a 
detailed analysis shows that there is an elongation of contours 
in the direction of the low-entropy tail. This is also confirmed 
by the ratio pressure map (Fig. |5]l which shows a clear devia- 
tion from radial symmetry in the SE sector and the presence of 
a secondary pressure peak. The entropy ratio map (Fig.|5]l con- 
firms the presence of low-entropy gas in the SE sector. These 
features in the ratio maps are still present even if we select the 
large scale center of symmetry as our center. 
It is apparent in our analysis that the regions where the ICM 
shows the most deviant features from the rest of the cluster are 
located in the SE sector 

In comparison to the typical scatter revealed in a representa- 
tive sample (Finoguenov et al. 2005), A3558 does not exhibit 
an anomalous amount of entropy or pressure substructures at 
large scales (r < 0.5rg()o). 



3.1. The cold front 

In order to characterize the edge seen in Chandra images, we 
have extracted the surface brightness profile across the edge 
(Fig. [8^): the discontinuity is apparent from the change of the 
slope at a projected distance ^ 100 arcsec from the cluster 
center 

We have extracted spectra from Chandra and XMM data in cir- 
cular annuli inside and outside the discontinuity. Fitting them 
with mekal models, we derive the temperature profile across 
the edge (Fig. [8]:): the gas is colder where the surface bright- 
ness is higher and therefore it is a cold front. Both instruments 
show the temperature jump across the discontinuity, but with 
the better statistics of the XMM-Newton observation, the tem- 
perature values measured by EPIC are more accurate. 
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Fig. 8. Upper left: Surface brightness profile obtained with Chandra across the discontinuity. Upper right: Electron density profile 
with the model described in Appendix A. Lower left: Projected temperature profiles across the cold front, obtained from ACIS 
spectra (red triangles) and EPIC spectra (black crosses). In order to reduce the effect of the EPIC PSF we did not extract the 
spectra in a slice of width 30 arcsec around the cold front. Lower right: Metal abundance profile across the cold front, using 
XMM-Newton spectra. 



Measuring the pressure jump across the c old front allows us 
to calcula te the speed of the moving cloud dLandau & Lifshitz 
1959t[viMilinin et alfcoOlh : 



y/(y-i) 



M < 1 



(> 

«M¥f (r-sr """'«> 1 



where Pst is the pressure at the stagnation point, that can be 
approximated by the pressure just inside the cold front. Pom is 
the pressure of the external, undisturbed gas, y is the adiabatic 
index and M is the Mach number of the motion of the cold 
cloud into the external gas. To calculate the pressure ratio, we 
derive the density jump from the fit of the surface brightness 
profile with the model described in Appendix A: we obtain 



= 1. 



0.13. We calculate the temperature jump 



with XMM-Newton since the statistical errors on the tem per- 
atures are smaller After deprojection dEttori et al.l 120021) . we 
obtain Tin/Tout = 0.85 ± 0.11. This leads to a pressure ratio 



Pst/ Pout - 1-60 ± 0.19, that corresponds to a subsonic motion, 
with a Mach number M - 0J8^q\1 
moving with a speed v = Mvs = 936^ 



with a Mach number M - 0.78^!J 1?. Therefore, the cloud is 

km s relative to the 
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gas outside the cold front. 

In the expression taken from iLandau & Lifshita (Il959h . Pout 
should be the pressure far upstream from the cold front, but 
unfortunately this does not apply easily to clusters since the 
pressure exhibits a strong gradient. For this reason, we have 
used for Pout the value in a region immediately outside the 
edge, which could be compressed and heated by the motion of 
the cold cloud. Therefore, we must consider the Mach number 
M - 0.8 a lower limit for the speed of the cloud. 
If the motion of the cold front were supersonic, we would 
expect the presence of a shock in front of the moving cloud. 
The shock would compress and heat the ICM leading to a 
second discontinuity in surface brightness and temperature. 
We have searched for the possible bow shock using the XMM 
observation, analyzing surface brightness and temperature 




300 300 400 600 600 



Fig. 7. Chandra image of the cluster in the band 0.8 - 3 k>^ v, 
adaptively smoothed using the CIAO task csmooth. The black 
curved line shows the position of the cold front, and the white 
elliptical sectors are used for the surface brightness profiles in 
Appendix A. Coordinates on the image are right ascension and 
declination. 



profile in elliptical annuli in the sector shown in Fig]?] Up 
to 200 arcsec from the edge of the cold front, there are no 
significant indications of a change in the slope in the surface 
brightness profile or of a temperature d iscontinuity. We can 
use Moeckel's method d VikhUnin et"al] (2001 ) to determine 
an upper limit to the Mach number of the moving cloud: a 
stand-off distance between the cloud and the shock of ^ 200 
arcsec corresponds to a Mach number M ^ 1.09. Therefore, 
we can conclude that the motion of the cloud is probably 
subsonic, with a Mach number in the range 0.78 - 1.09. The 
large indetermination on this result reflects the importance of 
the unavoidable assumptions on which it is based. Also the 
geometry and the symmetry of the system are very important 
and to prove this, we have used a simpler and less accurate 
model than that in Appendix A to fit the surface brightness 
profile across the edge, assuming spherical symmetry centered 
on the BCG. With this model we have found a lower density 
jump («/„/n„„f - 1.43 ± 0.04) and therefore a lower but 
consistent Mach number M - 0.53 + 0.07. 
The metalicity profile across the cold front obtained with 
XMM-Newton spectra (Fig. [8}J) could give us potentially 
important information on the origin of cold fronts. Large 
metalicity discontinuities are expected in some cold fronts 
due to merging where the core of an interacting subcluster 
is traveling through the ICM of a cluster poor in metals or 
through the outskirts of a cluster. The absence of a strong 
discontinuity in the metal abundance is consistent with several 



Fig. 9. Entropy map in the inner regions of the cluster The cir- 
cle radius is 113 arcsec, which corresponds to the distance be- 
tween the center and the cold front. The gas with lowest entropy 
is located west from the center, in the direction of the surface 
brightness discontinuity. Coordinates on the image are right as- 
cension and declination. 



possible scenarios: a "merging" cold front which travels 
through an ICM that h appens to have a simila r metalicity and 
also a "sloshing" core jMarkevitch et al ] |200lh . 
The mean metalicity inside the cold front, Z,„ - 
(0.41 + 0.03) Zq, is higher than the abundance in the gas 
outside, Zout = (0.31 ± 0.03) Z©, at a 3 cr level, but the profile 
is consistent with a smooth decreasing profile, like th ose 
observed in cool core clusters ( De Grandi & Molendil2001 ). or 
with a metalicity jump less than 0.1 (30-40 %). We conclude 
that the present data do not allow us to rule out any possible 
interpretation. 

The distribution of entropy inside the core of the cluster (inner 
120 arcsec) is intriguing: the gas with lowest entropy is located 
NW of the central AGN, near the discontinuity of the cold 



front (Fig. |9ll. In some numerical simulations (iHeinz et al 



2003; Ascasibar & Markevitch 20061) . the flow of gas past the 
moving cloud induces slow motions inside the cloud, which 
should transport gas from the center of the subcluster towards 
the interface. Our entropy distribution suggests that the low 
entropy gas located at the center of the cluster is moving 
towards the surface brightness discontinuity. 
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3.2. The anomalous SE sector 



As we have already outlined, there is a direction in A3558 
in which the ICM seems to have thermodynamic properties 
different from what we found in other sectors of the cluster: 
the SE regions. A low-entropy tail is observed in Fig. [TJl and 
Fig.[3}l, extending from the core toward the SE for ~ 500 kpc. 
Also the pressure is different in this sector, as seen in the ratio 
map (Fig.|5]): it is about a factor of two higher than predicted. 
The presence of unresolved clumps of denser matter could bias 
our estimate of pressure, since the density fluctuations enhance 
the cluster luminosity ove r what would be expe cted for a 
smooth, single-phase ICM dMathiesen et al. I ll999l) . From the 
test described in Appendix B, we conclude that the pressure 
excess is not an artifact due to unresolved clumps and that it 
must have a physical origin. 

Projected entropy and pressure map are perhaps not com- 
pletely reliable. To test the robustness of the detection of 
the anomalous thermodynamic properties of the SE sector 
we have compared these results to those obtained from an 
independent technique. In Fig. 13.21 we show the deprojected 
pressure and entropy profile in two opposite sectors of the 
cluster, obtained with spectral an alysis and the deprojection 
technique described in lEttori et alj ( 120021) . The entropy in the 
SE sector between 100 and 500 kpc is significantly lower than 
in the other direction, as predicted by the projected maps and 
the deprojected data are consistent with a pressure excess in 
these regions, even if at lower statistical significance. 
In order to study the metal distribution, we have extracted 
spectra from radial annuli in two sectors: one corresponding 
to the low entropy tail and the other to the rest of the cluster 
These profiles (Fig. fTTI) show that the abundance is higher in 
the regions of the low entropy gas than in the other regions. 
This suggests that the chemical history of the gas is different 
in the SE region with respect to other regions of the cluster. 
The anti-correlation between entropy and metal abundance 
is usually observed in the core of relaxed clusters. In the tail 
we find both low-entropy and high metal abundance (and also 
high pressure): this suggests that in the past this gas could 
have been at the center of a relaxed structure. Moreover, this 
also indicates that the processes that are responsible for the 
formation of this substructure did not significantly change the 
entropy of the ICM and therefore they are close to adiabatic. 
The results that come from the profiles are based on an 
underlying assumption: that the center from which we have 
extracted them (i.e. the peak of the surface brightness) is the 
center of the cluster We have some indications that the point 
we chose could be the center of the potential well of A3558 
(coincidence with the pressure peak ant the BCG), but we have 
already outlined that the center of the large scale symmetry 
does not coincide with this point (Sec. [3]). However the ratio 
maps with the new center still show the presence of high 
pressure and low entropy gas in these regions. Moreover the 
presence of gas with temperature similar to that of the core 
does not depend on the choice of the symmetry. 
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Fig. 10. Deprojected pressure (upper panel) and entropy (lower 
panel) profiles in the SE (green triangles) and NW (black 
crosses) sector, using XMM-Newton spectra. 

3.3. The core 

The metal abundance profile of the cluster (Fig. (TT] red boxes) 
shows a peak at the center, where also the surface brightness 
is peaked and the temperature drops (except than in the SE di- 
rection). These three features are always found together in cool 
core clusters: in A3558, we find a similar situation in a cluster 
which is not relaxed and which shows strong azimuthal vari- 
ations in the thermodynamic quantities, especially in the SE 
sector. However, all these features are not as apparent as in re- 
laxed objects. 

The deprojected central density is - (1.54±0.01)- lO^^cm"^ 
and the cooling time in the inner 40 kpc is tcooi - 5 ■ lO^yr, 
smaller than the Hubble time but larger than the typical values 
found in relaxed clusters dPeres et al..,1998.) 

4. Discussion 

Our analysis of the XMM-Newton and Chandra observations of 
A3558 has revealed many interesting characteristics. 
The cluster's dynamical history is probably more complicated 
than expected, since it cannot be classified easily as a merging 
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Fig. 11. Metal abundance profile in the low entropy tail (black 
diagonal crosses) and in the main cluster (red boxes) , using 
XMM-Newton spectra. 

or relaxed cluster. It has some features similar to those of cool 
core clusters such as: 

- An "intermediate" surface brightness peak at the center of 
the cluster (comparing Fig.|6]with the flux images of a sam- 
ple of nearby cluster we have found that A3558 surface 
brightness shows a clear peak at the center, while merg- 
ing clusters do not, even if this peak is less apparent than in 
cool core clusters); 

- A peak at the center in the m etal abundance profile, simila r 
to those in cool core clusters ( De Grandi & Molendil2001 ): 
in Figs.[8]l and[TT] we always find an increase of metalicity 
in the inner 60 kpc; 

- Presence of a BCG at the center and absence of other cD 
galaxies; 

- A small temperature drop (10 - 20%) at the center, except 
in the SE direction (Fig.[T]) where the temperature remains 
constant in the inner 200 kpc; 

- Cooling time smaller than the Hubble time in the inner 40 
kpc. 

However, it also has other properties that are more common in 
merging clusters: 

- Significant deviations from spherical symmetry in the ther- 
modynamic maps (as seen in Figs.[Tl[3]and also in the ratio 
maps in Fig. |5] the most deviating region is the SE sec- 
tor where we measure low temperature, low entropy, high 
pressure and high density); 

- Low entropy tail (with high metal abundance) in the SE 
sector (Figs. [U and O; 

- Pressure excess in the SE sector, detected with two inde- 
pendent methods (pressure ratio map, Fig.|5] and deprojec- 
tion. Fig. 13.21 1. 

Another puzzling featu re of this cluster is the absence of ex- 
tended radio emission dVenturi et alJl2000h . that would be ex- 
pected in massive merging clusters. 

These "strange" characteristics of A3558 make the interpreta- 
tion of its dynamical state a difficult task. From our data, we 



can exclude that the cluster is undergoing a major merger and 
that it is relaxed. 



4.1. Origin of the anomalous SE regions 

In Sec|3]we discussed that the SE sector of A3558 has "anoma- 
lous" thermodynamic properties with respect to the other sec- 
tors: low temperature (similar to that in the core), low entropy, 
high pressure and high metal abundance. Therefore, these re- 
gions are very important to study the dynamical state of the 
cluster In this section, we thus investigate different simple 
physical phenomena that could have lead to the formation of 
this anomalous structure. 



4.1 .1 . Projection effects 

One of the simplest scenarios to explain the SE sector 
"anomaly" is to invoke projection effects on the main cluster 
gas. In this hypothesis the low temperature measured in the SE 
sector could be due to a cold substructure superposed on the 
line of sight of the cluster, while A3558 would have spheri- 
cally symmetric distributions of the thermodynamic quantities 
and would look "relaxed". 

To test this hypothesis we fitted the spectra extracted in a cir- 
cle in the "tail" (Fig.|5]l with a 2-temperature model but the fit 
does not improve significantly and the temperature of the sec- 
ond component is very similar to that of the first. These results 
show that the proposed scenario cannot explain the observed 
data: even if there were a second thermal component due to a 
projected substructure, the temperature of the ICM of the main 
cluster would remain lower than in the other sectors and we 
would still need another mechanism to explain this. 
Since the anomaly is located in the infall region, where the fil- 
ament connects A3558 to the other members of the core of the 
Shapley supercluster, a possible superposed substructure could 
be the filament itself. Therefore, we fitted the spectra with two 
thermal components, fixing the value of the second tempera- 
ture in the range 1 - 1.5 keV (Kull & Bohringer 1999). Again 
the quality of the fit does not improve with this model and from 
the ratio of the normalizations of the components, we find that 
a filament can contribute less than 1% to the emission of the 
region. 

4.1.2. An embedded subcluster 

Since the pressure distribution is a tracer of the shape of the 
potential well of a cluster, the pressure excess in the SE sector 
may indicate the presence of the dark-matter halo of a small 
structure which is interacting with the main cluster. This sce- 
nario is consistent with the low temperature and entropy and 
the high metalicity measured in this region. 
We have characterized the subcluster, using a model similar 
to the one described in Appendix B, assuming that region 2 
in Fig. [14] is a substructure completely filled by high density 
gas. Fitting the spectrum with a two-temperature model (one 
for the substructure and the other for the emission along the 
line of sight), we find that the substructure has temperature 
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kT = (3.82+0.07) keV anddensityn = (2.54+0.03)-10-^cm-3 
Q, from which we can derive its gas mass Mga^ - 3.89 ■ 10'' 
Mq. Assuming hydrostatic equilibrium, we have calculated the 
total mass contained in the pressure peak: it is in the range 
(3.5 - 7) ■ 10'2 Mo, lower by a factor of 3-5 than the mass con- 
tained in the same radius at the center of the cluster. Since this 
model is based on strong assumptions, especially the choice 
of the surface brightness peak as the center of the cluster, we 
should treat the parameters that we have derived as rough esti- 
mates. 

A mass ~ 5 ■ 1O'^M0 is comparable with the mass of a large cD 
galaxy or of a small group, which usually have a virial temper- 
ature lower than the value that we derived for the subcluster 
This discrepancy can be explained since we are not observ- 
ing an isolated group, but we are describing the properties of 
a subcluster embedded in a massive cluster. The subcluster has 
probably entered the main cluster from SE, from the "filament" 
connecting A3558 to the SC group. This is the preferential di- 
rection of merging and the SE sector is therefore the infall re- 
gion of matter in the potential well of the cluster 
There are two possible critical points of this model: the absence 
of correlation between the pressure excess and optical data and 
the lack of heating and visible hydrodynamic effects associated 
with the infall of the subcluster 



Com paring our pressure map and optical data (IBardelli et al 



19981) . we find no clear indication of a particular concentration 



of galaxies corresponding to the pressure excess. However this 
is a weak objection since we do not expect this substructure to 
be associated with a large number of galaxies, because of its 
low mass. 

The other objection is more serious: we would expect a merg- 
ing subcluster to produce disturbances visible in the X-ray im- 
age and in the thermodynamic maps, while there is no evidence 
of heating in the infall region, which is furthermore the coldes t 
sector of the cluster Recently, Ascasibar & MarkevitchI ( 2006 ) 
have shown a set of simulations where the merging between 
a massive cluster and even small subclusters produces global 
disturbances in the ICM visible in the X-ray image long after 
the first core passage. However , in another set of simulations, 
Ascasibar & Markevitch (l2006l) show that if the infalling sub- 



cluster has no gas there are no visible disturbances in the X-ray 
image, except the sloshing of the core (that we will discuss 
in Sec l4.2b . This could explain the lack of heating in the in- 
fall region of A3558: if the gas of an infalling subcluster has 
been completely stripped moving through the dense environ- 
ment of the filament before reaching the cluster, we should 
treat the subcluster as "dark-matter" only and therefore we 
should not expect large hydrodynamic effects. However in this 
scenario, the subcluster would not contribute its low tempera- 
ture and entropy ICM in the SE regions of the cluster. It can 
only perturb the gravitational potential and induce some hy- 
drodynamic eff'ects that could move the low-entropy gas from 
the core of the main cluster to the outer regions (in the case 
of A3558 up to 500 kpc from the surface brightness peak) 
d Ascasibar & Markevitchll2006t) . 



cal. 



; errors reported for temperature and entropy are only statisti- 



4.1 .3. Stripping due to tine motion of a cold cloud 

From the analysis of the cold front (Sec. 13.1b . we know that 
A3558 hosts a cloud of dense and cold gas which is currently 
moving NW. During its motion through the cluster it could have 
lost part of its high density and low temperature gas ("core- 
like" ICM) because of the stripping from the surrounding envi- 
ronment. Indeed, the tail-shape of the low entropy region (Fig. 
[U) suggests that some stripping has occurred. 
This mechanism can be effective only if the ram pressure is 
large enough to overcome the internal pressure of the subclus- 
ter, i.e. in the condition 



Penv^ Pc Penvn 



(3) 



where peni- is the density of the environment in which the cloud 
is moving with velocity v, Pc is the internal pressure of the 
cloud and f e„,. the pressure of the surrounding environment. 
We have calculated the velocity necessary to satisfy Eq. |3]for a 
cold cloud that travels through the ICM of A3558: 



y(r) = 400km s^'keV 



-1/2 



n^Tc - n(r)T(r) 
n(r) 



1/2 



(4) 



where n(r) and T(r) are the properties of the ICM, as a function 
of the distance from the cluster center. In the hypothesis that 
the moving cloud can be described by the properties of the 
ICM inside the cold front, we derive that the required speed is 
consistent with the value obtained from the analysis of the cold 
front only at distances smaller than 150 kpc from the cluster 
center. The low entropy gas located at the end of the tail (500 
kpc) cannot have been stripped and deposited "in situ" during 
the motion of the cloud through the cluster, since the stripping 
can be effective only in the inner regions where the density of 
the environment is large. 



4.1 .4. A3558 former core 

One interesting feature of A3558 is the fact that the center of 
the large scale isophotes does not coincide with the surface 
brightness and pressure main peak but is located ~ 70 arcsec 
SE of the BCG. This effect has been observed in many clusters, 
but it is interesting that in A3558 the X-ray centroid is located 
in the pressure excess observed in the ratio map (Fig.[T2]l. 
This coincidence may have an interesting physical meaning: 
the large scale isophotal centroid could indicate the original 
position of the core before the processes that have lead to 
the formation of the cold front (Sec 14.2b . We can assume 
that the center of the cluster was located at the actual po- 
sition of the X-ray large-scale centroid, then that the core 
started to oscillate producing a cold fr ont ( Tittley & Henriksenl 
I 2OO5I: lAscasibar & MarkevitchI l200d) and therefore distorting 
the isophotes in the inner regions, while the external part of 
the cluster is not perturbed. The low entropy, high pressure and 
high metal abundance gas located near the former center would 
be part of the ICM of the core. 

The presence of low-temperature and -entropy gas SE of the 
surface brightness center (Fig. [TST i can be explained only 
taking into account hydrodynamic processes. Figure 7 in 
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Ascasibar & MarkevitchI (120061) shows a similar situation even 
if at smaller scales: the ICM of the sloshing core is dragged 
outward by the motion of the surrounding gas. 

4.1 .5. Summary 

We have proposed several simple scenarios to explain the 
anomalous properties of the SE sector Projection effects are 
ruled out by our analysis, the presence of an embedded sub- 
cluster requires very strong assumptions on the nature of the 
substructure and the stripping deposition in situ by an infalling 
subcluster would require higher velocities or densities. It is ap- 
parent from our analysis that all simple scenarios fail to ex- 
plain the properties of the SE region s and that hydrodynamic 
proces ses, such as those described in lAscasibar & Markevitch 
(l2006h . should be taken into account. 



4.2. Origin of the cold front 



In Sec. 13.11 we have shown that A3558 hosts a cold front, 
moving transonically with a Mach number M ^ 0.8. The 
metal abundance profile across the brightness discontinuity is 
consistent with a smooth decreasing profile and with a small 
(30 -40%) jump. 

Several mechanism have been proposed in the literature for 
the origin of cold fronts: they could b e the remnant of a 
merging subcluster, as in A3667 (Vik hlinin et al.l 1200 lb . or 
the core o f the cluster "sloshing" in the potential well, as 



in A 1795 (M arkevitc h et al.l 12001 
Ascasibar & M arkevitc hll2006h~ 



Tittlev & Henriksen 2005 



If the cold front that we have detected is the remnant of a sub- 
cluster that is falling on A3558 with a subsonic motion, we 
would expect to recognize also the parent cluster's core, and 
this is not clearly detected either in the X-rays or in the op- 
tical. To support this hypothesis we have to make one of two 
strong assumptions: either that the moving subcluster is super- 
posed on the line of sight to the core of the main cluster or that 
the main cluster core has been almost completely destroyed by 
the merging subcluster. Moreover, since we do not observe a 
metalicity discontinuity across the cold front, we also have to 
assume that the metal abundance distribution of the subcluster 
must be similar to that of the main cluster 
The other possible origin of the cold front is the "sloshing" 
scenario, which can describe two phen omena: the oscillation 
of the gas component of the core (M arkevitch et al.l 1200 lb 
and the oscillation of the whole core, dark matter and gas 
in response to the off-a xis merger between t he clu ster and a 
group, as suggested by 



Ascasibar & Markevitch 



ttlev & HenriksenI (l2005b . Recendy 
( 2006h lave shown that these two 



cases are not separate, since a perturbation of the gravitational 
potential induces motion of the gas, which quickly decouples 
from its dark matter compone nt. However, we c o mpare our re- 
sults also with the models of Markevitch et al.l (120011) and of 
Tittlev & Henrikse"nl (2005), because they predict simpler ob- 
servational properties, which can be more easily compared with 
real data. 



from its dark matter component and oscillates in the potential 
well of the cluster. Since in A3558 the pressure peak of the 
cluster (i.e. the bottom of the potential well) coincides with the 
peak of the surface brightness (i.e. the center of the gas cloud), 
as we can see in Fig. [T] we conclude that, in this scenario, the 
moving gas cloud is now in the center of its oscillation and 
therefore in the point of maximum velocity. This would also 
explain the high Ma ch number (M ^ 0.8), compared to other 
sloshing cold fronts (iMarkevitch et al ] l200lb . 
If the sloshing is due to the oscill ation of the gravitational 
potential dTittlev & Henriksenll2005 ). the coincidence between 
pressure and surface brightness is expected and we cannot de- 
termine the position of the oscillating b ody with respect to 
the ce nter of its oscillation. In their paper. iTittlev & Henriksen 
(l2005l) lighlight three main observational signatures of this 
kind of cold front, but we could verify only the third (the gra- 
dients of compressed density must point through the center of 
the cluster) because we detected only one cold front. 
The environment of superclusters is ideal to observe the slosh- 
ing class of cold fronts, since the merging rate is higher than 
in field clusters. We have identified two groups that may be re- 
sponsible for the perturbations in the potential well: the "sub- 
cluster" (if we interpret the anomalous SE sector as a merging 
substructure as in Sec. 14.1.2b and the group SC 1327-312. In 
either case, to be responsible for the oscillation of the core, 
these structures must have already passed the point of closest 
encounter with A3558. It is difficult to suppose that the low- 
mass subcluster has already reached the point of closest en- 
counter with the core of A3558: we have to assume that it has 
entered with a non-zero angle with respect to the plane of the 
sky and it is now behind or in front of the cluster. In this sce- 
nario, the core would have already completed half of its oscil- 
lation. A past closest approach between A3558 and SC 1327- 
312 seems more plausible. The group could have passed north 
of the cluste r to reach its a ctual position. A similar scenario is 
proposed bv lBardelli et al. I ([1998): they suggest that A3558 and 
another cluster (whose remnants are A3562, SC 1329-313 and 
SC 1327-312) colUded in the past. 

It is difficult to compare our data more quantita tively with the 
model proposed by Tittl ev & Henrikse"nr ( 2005 ). because their 
results are based only on two sets of simulations with fixed 
mass ratios (30:1 and 10:1) and none of our candidate groups 
has a similar mass ratio with respect to A3 558 (~ 100 : 1 fo r 
the subcluster and ~ 4 : 1 for SC 1327-31 2, [Ettori et al. 1997b. 



In the model of Markevitch et al.l (120011) . the gas has decoupled 



Moreover, as outlined by l Ascasibar & Markevitch (2006), 
complex hydrodynamic processes should be taken into ac- 
count. 

Comparing our results with the possible cold front models does 
not allow us to rule out completely any possible interpretation. 
However we prefer the "sloshing" model, since the interpreta- 
tion of the cold front of A3558 as the remnant of a merging 
subcluster would require too much coincidence. 



4.3. The global picture 

In this section we attempt to develop a unified picture and to in- 
terpret the dynamical state of the cluster. We have excluded that 




Fig. 12. The crosses show the position of the center of large scale symmetry. Upper left: EPIC flux image with the elUpses used 
to approximate contours. Upper right: temperature map , Lower left: pressure ratio map and Lower right: pseudo-entropy map. 
The green cross shows the position of the centroid of X-ray large scale isophotes. 



the cluster is relaxed and we have shown that all the main fea- 
tures of the cluster are consequences of interaction. However, 
we still have to discriminate between two classes of merging 
scenarios: one where the observed features are a direct effect 
of the interaction and the other where we see the consequence 
of hydrodynamic processes that follow the perturbation of the 



gravitational potential. 

A possible scenario is the infall of a subcluster, that has en- 
tered A3558 from SE and is currently moving towards NW. 
This merging will produce, as direct effects, the cold front (sep- 
aration edge between the ICM of the subcluster and that of 
the main cluster) and the stripped tail. However, we calculated 
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(Sec. |4.1.3t that a higher velocity would be required for the 
stripping to form the tail. Moreover, a subcluster falling in the 
potential well of A3558 would have acquired a kinetic energy 
much greater than the value that we derive from the velocity of 
the cold front. This "missing energy" should have been trans- 
formed into thermal energy, but we have no evidence of heat- 
ing. This merger should have produced disturbances visible in 
the morphology of the cluster, while the presence of a large 
scale symmetry disfavors this scenario. We have already shown 
(Sec. 14.2b that the interpretation of the cold front as a merging 
subcluster would require too much coincidence: the superposi- 
tion between the subcluster and the core and similar metalicity. 
The main features of A3558 can also be interpreted as an in- 
direct effect of an interaction. The "sloshing scenario" can ex- 
pl ain not only the cold front, but also the tail. As explained 



Ascasibar & Markevitch (l2006l) the cool gas would be dis- 



placed from the potential minimum and, because of the ram 
pressure of the surrounding gas, part of this cool gas would 
expand adiabatically and move outward. This scenario is con- 
sistent with the observed large scale symmetry, whose center 
does not coincide with the surface brightness peak. Finally the 
absence of current strong interactions with other structures can 
explain the lack of extended radio emission and the presence of 
a sort of cool core. 

Our analysis favors a scenario where the gravitational potential 
of A3558 has recently been perturbed by a merger that has in- 
duced oscillations of the core without destroying the structure 
of the cluster A qualitative comparison between our data and 
the simulatio ns shows that hydrodyna mic effects, such as those 
described in Ascasibar & MarkevitchI (i2006i) . may indeed play 
an important role in the physics of the ICM. 




Z 



5. Conclusions 

Our analysis has shown that A3558 cannot be considered either 
a relaxed cluster, even if it has a sort of "cool core" in its inner 
regions, or a merging one. It hosts a cold front, moving NW, 
which is probably due to the oscillation of the gravitational 
potential. During these oscillations, hydrodynamic processes 
have produced a "tail" of low-entropy, high-pressure and metal 
rich ICM. The lack of current strong interactions with other 
structures can explain the absence of extended radio emission. 
However, something has perturbed the dynamical state of 
A3558. This "perturbing factor" is probably now outside 
A3558: a past merger with SC 1327-312 or with a more 
massive cluster whose remnants are A3562, SC 132 9-313 
and SC 1327-312 (as suggested bv lBardelli etal]|l998b . may 



have perturbed the gravitational potential. Indeed, in the 
processes related to th e formation of a "supercluster" structure 
(lAkimoto et al.ll2003l) . there are frequent interactions that can 
have induced the "sloshing" of the core. 



Fig. 13. Geometrical model used to derive the density discon- 
tinuity across the cold front. C is the centroid of the surface 
brightness large scale contours, SB is the surface brightness 
peak, D is the distance between the centers, b is the projected 
distance, Z is the position along the z axis (perpendicular to the 
plane of the sky), and rc are the distances calculated from the 
surface brightness peak and the centroid, respectively. 

Appendix A: Fit of tlie SB profile across the cold 
front 

The density jump across the cold front can be derived from the 
surface brightness profile across the discontinuity. We assume 
the density profile to be composed of two [5 models 



«i ^1 + (t;) j outside 



inside the cold front 



(5) 



and we derive the parameters of this model (n^, ri, a, rih, r2, P) 
and the projected distance of the edge xq (calculated from the 
brightness peak) from the fit of the surface brightness profile, 
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which can be expressed as the integral of the emissivity along 
the line of sight: 



1/2 p 

Jo 



rr{z)dz, 



(6) 



where K is a constant. Since the dependence of the surface 
brightness on the temperature is weak, we do not consider its 
variation along the line of sight and we can leave it outside the 
integral. 

The choice of the symmetry is very important: the profile inside 
the cold front should be centered on the center of curvature, 
which in our case coincides with the surface brightness peak, 
while the outer/? model should be centered on its own center of 
symmetry. From the analysis of the XMM-Newton observation 
(Sec. |3]l, we know that the large scale symmetry is elliptical 
and the center does not coincide with the surface brightness 
peak but is located 70 arcsec SE. This geometry is reproduced 
in Fig. [13] where z is the line of sight, and rc are the distances 
calculated from the surface brightness peak and the large scale 
symmetry center respectively and D is the distance between the 
centers, that we assume equal to the projected distance 70 arc- 
sec. 

To match the symmetry of the outer region, we have extracted 
the surface brightness profiles in four small elliptical sectors 
shown in Fig.|2l in each of these sectors we can locally assume 
that the symmetry is spherical, and we can write the surface 
brightness as: 



2 = 2B 



Jo 



dz 



(7) 



where is the distance calculated from the large scale center 
and B = KT^^V^. Changing variables (/ = z^/(rl + b^)), this 
expression can be written as: 



/ I Jo 



dy 



(1 +3;2)3/^' 



(8) 



where be = b + D is the projected distance from the large scale 
symmetry center. 

Since we are in the conditions b^ » r2 the expression in 
Eq. [8] can be approximated by a power law S = Bb^^, where 
7 = 6y3 - 1. We fit the profiles and we derive j - 1.36 + 0.09 
and the normalization. 

The expression for the surface brightness inside the edge is 
more complicated since we must consider two components: 
one for the emissivity of the inner region integrated between 

and -^Xg - b^ and the other for the emissivity of the outer 
region between - b^ and oo: 



2A 



dz 



+ 2B 



r dz 



(9) 



where A = Q.^2KT^'\l and r, and are the distances cal- 
culated from the surface brightness peak and the large scale 
centroid respectively and xq is the curvature radius of the cold 
front. Under the same hypothesis that we have used in the pre- 
vious case (b » r2), the second part of Eq.|9]can be expressed 



B 



f 



dy 



(1+/ 



Jo 



dy 
(1+y^y'H 



(b + D)' 



,(-6/S+I) 



(10) 



where y = - b^/(b + D) and b is the projected distance 
from the surface brightness peak. 

The first term in Eq. |9] cannot be approximated with a power 
law because we are not in the case b » r\ : 



2Ari 



2\ 



-3Q-+1 



dy 



(i+y^y 



2\io ■ 



Fitting the surface brightness profile in the inner region, ex- 
tracted in the spherical sector shown in Fig. |71 with the func- 
tion S = 2i H- ^2, leaving the parameters of S2 fixed to the best 
fit value for the outer region, we derive the parameters n^, ri 
and a. Evaluating the y6 model at r - 110 we derive the den- 
sity n,„ at the position of the discontinuity and the density jump 
nin/noui - (1.43 ± 0.04). The best fit model for the density is 
shown in Fig.[8j). 

Appendix B: The clumping test 

To test if the observed excess could be due to clumping we have 
extracted the spectra in the region where the surface brightness 
excess is greatest (the white circle in Fig.|5]l, and fitted it with 
a two temperature model. To account for the projection effects 
we had to assume some geometrical and physical conditions. 
We consider two identical cylinders (Fig. fl4l i of radius R and 
height h - 2R, with the bases perpendicular to the line of sight 
(the bases corresponds to the circles in Fig.|5] the white circle is 
the base of cylinder 2 while the black one is that of cylinder 1). 
We assume that the cluster is spherically symmetric around the 
surface brightness peak , with the only dishomogenity located 
in cylinder 2 and filling its volume with a filling factor (p. The 
emission in direction 1 is due to the homogeneous gas with 
density ni of cylinder 1 and to the gas on the line of sight with 
density rieir): the XSPEC normalization of the mekal model 
that we have used to fit the spectrum can be expressed as 



/-*oo 

^2k I I 

Jh/2 



nl{r)Adr + kn\V, 



where A is the area of the base, V 
cylinder, k is the XSPEC constant 

10-14 



(11) 



Ah is the volume of the 



47r(D„(l+z)2)' 



(12) 



and Da is the angular size distance to the source (cm). The value 
of the gas density alb ^ 200 kpc, ni = 7.52 ■ lO^'^cm--', has 
been derived from the deprojection in spherical sectors. 
The emission of cylinder 2 has two thermal components: one 
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direction 2 direction 1 



observer 

Fig. 14. Geometrical model adopted for the clumping test.R is 

the radius of the circles from which we extract the spectra, h is 
the width of the region in which we assume there are clumps 
(along the line of sight), b and r are the projected and three- 
dimensional distances, respectively, calculated from the surface 
brightness peak. 

due to the clump with density ric and the other due to the ho- 
mogeneous gas with density rii, filling the remaining volume. 
The normaUzation due to this second component (^2) can be 
expressed in terms of the normalization in direction 1, A'^i : 

N2 = 2k[ nlir)Adr + knjil - 4>)V = Ni - knl(pV. (13) 

Jh/2 

We have fitted the spectrum from direction 2, with a two- 
temperature model varying the temperature of the "clump" 
component and its filling factor. The best fit combination 
of these parameters, Tc = 3 keV and (f> = 0.1, leads to a 
density in the clump which is about ten times higher than the 
density of the other component. Since the temperature ratio 
is much smaller (Tu/Tc = 1.97), the clump component is not 



in pressure equiUbrium with the surrounding medium. Even 
with greater filling factors the pressure equilibrium can never 
be achieved, unless we force the temperature to unrealistically 
low values which worsen significantly the quality of the fit. 
This analysis indicates that the pressure enhancement observed 
in the ratio maps is not an artifact due to unresolved clumps 
and that it must have a physical origin. 
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